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Abstract 

Cohomology and cohomology ring of three-dimensional (3D) objects are topo- 
logical invariants that characterize holes and their relations. Cohomology 
ring has been traditionally computed on simplicial complexes. Nevertheless, 
cubical complexes deal directly with the voxels in 3D images, no additional 
triangulation is necessary, facilitating efficient algorithms for the computa- 
tion of topological invariants in the image context. In this paper, we present 
formulas to directly compute the cohomology ring of 3D cubical complexes 
without making use of any additional triangulation. Starting from a cubical 
complex Q that represents a 3D binary-valued digital picture whose fore- 
ground has one connected component, we compute first the cohomological 
information on the boundary of the object, dQ by an incremental technique; 
then, using a face reduction algorithm, we compute it on the whole object; 
finally, applying the mentioned formulas, the cohomology ring is computed 
from such information. 

Keywords: Cohomology ring; cubical complexes; 3D digital images. 



1. Introduction 

Many computer application areas involve topological methods which usu- 
ally mean a significant reduction in the amount of data. Homology is an al- 
gorithmically computable topological invariant that characterizes an object 
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by its "holes" (in any dimension). Informally, holes of a 3D-object are its 
connected components in dim. 0, its tunnels in dim. 1 and its cavities in dim. 
2. Cohomology is a topological invariant obtained by an algebraic duality 
of the notion of homology. Although the formal definition of cohomology is 
motivated primarily by algebraic considerations, homology and cohomology 
of 3D objects are isomorphic, that is, they provide the same topological infor- 
mation. Nevertheless, cohomology has an additional ring structure provided 
by the cup product (denoted by ^). The cup product can be seen as the way 
the holes obtained in homology are related to each other. For example, think 
of the torus, and the wedge sum of two loops and a 2-sphere. Both objects 
have two tunnels and one cavity; but the cavity (7) of the first object can be 
decomposed in the product of the two tunnels (a and j3), that is, a ^ (3 = 7; 
the cavity of the second object cannot (see Fig. [1]). This information would 
contribute to a better understanding of the degree of topological complexity 
of the analyzed digital object, and would shed light on its geometric features. 

In 0, IH, a method for computing the cohomology ring of 3D binary 
digital images is stated. In those works, the cohomology ring computation 
is performed over the (unique) simplicial complex associated to the digital 
binary- valued picture using the 14- adjacency. However, one could assert 
that a more natural combinatorial structure when dealing with 3D digital 
images is the one provided by cubical complexes. One way to compute the 
cohomology ring of a cubical complex Q is to convert it into a simplicial 
complex Kq by subdividing each cell and applying the known formulas for 
computing the cohomology ring of Kq. In this paper, we present formulas 
to directly compute the cohomology ring of 3D cubical complexes without 
making use of additional triangulations. Besides, we describe a strategy 
to tackle the cohomology ring computation on a 3D binary-valued digital 
picture. This paper extends a preliminary version (see j^j). 

The paper is organized as follows. In Section [21 we recall the concept of 
AT-model and extend it to general polyhedral cell complexes; given the AT- 
model for a polyhedral cell complex, we give the formulas of a new AT-model 
obtained after a subdivision; this result will be the key to prove the validity 
of the formulas of the cohomology ring of cubical complexes. In Section [21 
formulas for computing the cohomology ring of 3D cubical complexes are 
established. Section H] is devoted to describe the process to obtain an AT- 
model of a 3D digital image that provides the ingredients for the cohomology 
ring computation. Finally, some conclusions and plans for future are drawn 
in Section O 
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Figure 1: On the left, a hollow torus and its two tunnels. On the right, the wedge sum of 
a 2-sphere and two loops. 



2. AT-models for Polyhedral Cell Complexes 

In this section, we recall first the concept of AT-model (see 0,01) for a cell 
complex, which consists of an algebraic set of data that provides homological 
information. Given the AT-model for a polyhedral cell complex, we show the 
formulas of a new AT-model obtained after a subdivision; this result will be 
the key for the formulas of the cohomology ring of cubical complexes in the 
next section. 

Since we are working with objects embedded in M 3 , the homology groups 
are torsion-free 0, ch.10], so computing homology over a field is enough to 
characterize shapes [HI, p. 332]. This fact, together with the isomorphism 



(over any field) between the homology and cohomology groups [ll|, p. 320], 
enables us to consider Z/2 as the ground ring throughout the paper. 

A polyhedral cell complex P in I 3 , is given by a finite collection of cells 
which are convex polytopes (vertices, edges, polygons and polyhedra), to- 
gether with all their faces and such that the intersection between two of 
them is either empty or a face of each of them. A proper face of a G P is a 
face of a whose dimension is strictly less than the one of o. A facet of a is 
a proper face of a of maximal dimension. A maximal cell of P is a cell of P 
which is not a proper face of any other cell of P. Observe that if the cells of 
P are n-simplices, P is a simplicial complex (see [IH); in the case that the 
cells of P are n- cubes, then P is a cubical complex (see 0]). A g-cell of either 
a simplicial complex or a cubical complex can be denoted by the list of its 
vertices. 

For any graded set S = {S q } q , one can consider formal sums of elements of 
S q , which are called q- chains, and which form abelian groups with respect to 
the component- wise addition (mod 2). These groups are called q-chain groups 
and denoted by C q (S). The collection of all the chain groups associated to 
S is denoted by C(S) = {C q (S)} q and called also chain group, for simplicity. 
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Let {s\, . . . , s m } be the elements of S q for a fixed q. Given two g-chains 
c i = Yn=i a i s i and c 2 = YhLi where G Z/2 for i = l,...,m, the 

expression (ci,c 2 ) refers to YlT=i a i ' A e Z/2. For example, fixed z and j, 
the expression (ci, Sj) is a« and (sj, Sj) is 1 if i = j and otherwise. 

The polyhedral chain complex associated to the polyhedral cell complex 
P is the collection C{P) = {C q (P),d q } q where: 

(a) each C q (P) is the corresponding chain group generated by the 5- cells 
of P; 

(b) the boundary operator d q : C q (P) — > C q -i(P) connects two immediate 
dimensions. The boundary of a g-cell is the formal sum of all its facets. 
It is extended to g-chains by linearity. 

For example, consider a triangle (vi,Vj,Vk) with vertices Vi, Vj, Vk- The 
boundary of the triangle is the formal sum of its edges, that is, d 2 (vi, Vj,Vk) = 

(Vi,Vj) + {Vj,V k ) + (Vi,V k ). 

Given a polyhedral cell complex P, an algebraic-topological model (AT- 
model jl, III) for P is a set of data (P, H, f, g, 0), where H is a graded subset 
of P and /, g, are three families of maps {f q : C q (P) — > C q (H)} q , {g q : 
Cg(H) -t C q (P)} g and {(f) q : C q (P) -> C q+ i(P)} q , such that, for each q: 

(!) fq9q = idc q (H), (pq-id q + d q+1 (j) q = id Cq ( P ) + g q f q , f q -\d q = 0, d q g q = 0; 

(2) <p q+l 4> q = 0, f q+l 4> q = 0, cf) q g q = 0. 

As a result, the chain group C(H) is isomorphic to the homology (and to the 
cohomology) of P. In particular, the number of vertices of H coincides with 
the number of connected components of P, the number of edges of H with 
the number of tunnels of P and the number of 2-cells of H with the number 
of cavities of P. Fixed q, for each a G H q , g q (cr) is a representative cycle of 
a homology generator of dim. q. Define an homomorphism a*f q : C q (P) — > 
Z/2 such that if /i is a g-cell of P then 

^fqdA ■= (o-, fq(fi)) mod 2. 

Then, cr*f q is a representative cocycle of a cohomology generator of dim. q. 
An isomorphism between C(H) and the homology (resp. cohomology) of P 
maps each a G H to the homology class represented by ^(cr) (resp. the 
cohomology class represented by cr*f q ) (see [4,0]). 

From now on, we will omit subscripts on behalf of simplicity. 
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Figure 2: a) An abstract cubical representation Q of the hollow torus; b) the paths "/( VitVQ ) ; 
c) the "path" c e , for any edge e of Q. 



Example 2.1. Let Q be an abstract cubical representation of the hollow 
torus given in Fig. An AT -model for Q is the set of data (Q, H, f, g,<p) 
given in the following table: 



Q 


f 


<\> 


H 


9 


vo 


v 





vo 


vo 


Vi, i = 1, . . . ,8 


vo 


7(vi,Uo) 






a i; i = 1,2 


bi 








h, t = l,2 


bi 





k 


Oti 


any edge b ^ Oi,bi 











c 


c 





c 


p 


any square a ^ c 













where a x e {(t> 3 , fe), K, «b)}; «2 e {(t>i, t> 7 ), (t> 2 , t> 8 )}/ &i = (v ,v 2 ); b 2 = 
(v , Va); c = (v , v 2 , W4, v 8 ); 7(t, ijVo ) £/ie on/y pa£/i m Q /rom t> j to t>o 

f/or example, j( V7jVo ) = (v 5 ,v 7 ) + (v u v b ) + (v ,fi)j; given an edge e o/ 
Q, c e is the sum of the squares that correspond to the path starting from e 
and following the arrows in Fig. [He (for example, C(„ ljt , 7 ) = (v ,Vi,V4,v 7 ) + 
(^3,^5,^4,^7) + (v ,Vi,v 3 ,v 5 )); representative cycles of homology generators 
are the vertex v , the tunnels a.\ = (v ,vi) + (vi,v 2 ) + (v ,v 2 ) and a 2 = 
(vq, U3) + (vs, V4) + (vq, U4) and i/ie cavity /3 which is the sum of the 9 squares 
ofQ. 

An algorithm for computing AT-models for polyhedral cell complexes ap- 
pears for example in 0, In fact, in those papers, the algorithm is designed 
for simplicial complexes but the adaptation to polyhedral cell complexes is 
straightforward. That algorithm runs in time 0(m 3 ) where m is the number 
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Figure 3: A subdivision of a square a in two triangles ai and a2- 



of cells of the given polyhedral cell complex. If an AT-model (P,H,f,g,(f>) 
for a polyhedral cell complex P is computed using that algorithm then it is 
also satisfied that if a G H then f(a) = a and a G g(a). 

Let P be a polyhedral cell complex. Fixed q, we say that a g-cell a G P 
is subdivided into two new g-cells a\ ^ P and «2 ^ P by a new (q — l)-cell 
e ^ P (see, for example, Fig. [3]) if: 

(a) e is a facet of «i and a 2 ; 

(b) «i U «2 = «; 

(c) CKi fl a 2 = e. 

The following lemma establishes how to obtain a new AT-model for P 
after subdividing by a (g — l)-cell e a g-cell a of P into two new g-cells «i 
and « 2 . This result will be used to prove the equivalence between the new 
formula to compute the cup product on cubical complexes and the classical 
one on simplicial complexes. 

Lemma 2.2. Let (P, H, f, g, <p) be an AT-model for a polyhedral cell complex 
P computed using the algorithm given in 0/ Let a be a q-cell, which is 
subdivided into two q- cells a.\ and a 2 by a new (q — I) -cell e. Let H' := (H \ 
{a}) U {ai} if a £ H, and H' := H otherwise; P' := (P \ {a}) U a 2 , e}; 
and & the boundary operator of P' given by: d'(c) := <9(c) + (a,d(c))(a + 
a.\ + a 2 ) for any c G P'\{e, ai, a 2 }. Denote <9(«i) + e by A, and <9(a 2 ) + e by 
B. Then, the set (P',H',f',g',(j)') is an AT-model for P' , where f , g' and 
0' are given by: 

. := f(a) + (aj(a))(a + a 1 ); f\a 2 ) := 0; /'(e) := f(A)(= f(B)); 

f'W) '■= f( a ) + («>/( cr ))(« + oti), for any a G P' \ {a 1 ,a 2 ,e}; 

• 0'(ai) := 0(a); <p'(a 2 ) := 0; 0'(e) := a 2 + 0(P) + (ct, 0(P))(a + ai + a 2 ); 
0'(<r) := 0(a) + (a, <j)(a))(a + cm + a 2 ), for any a G P' \ {cm, a 2 , e}; 
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• If a G H, g'{a\) : = g(a) + a + a.\ + a 2 ; 

g\l) : = din) + («,5'(7))(« + «i + ot 2 ), for any 7 G H' \ {«i}. 

ProofWe have to check that (P', H', f, g', 0') is an AT-model for P' . We 
will only check that f'g' = id and id + g'f = 0'9' + & <f> . The rest of the 
conditions are left to the reader. 

Let 7 G H', 7 ^ cki, and <r G P'\{qi, a 2 , e}. 

• /V = id: 

/V(7) = f (^(7) + («^(7))«) + («^(7))f(«i) = fg(l) + (a,g( 1 ))f(a) + 
(a,g(-y))(a,f(a))(a + ax) + (a, g(-y))(f(a) + (a, f(a))(a + an)) = 7. 

If a G B; /V(ai) = f (»(a)+a)+/' (ai) = f(g(a)+a)+f(a)+a+cn = 

• jd + ^'/' = + 9'0': 

ax + g'f'{oL\) = oil + g'(f(a) + (a, f(a))a) + (a, f(a))g'(ai) = an + 
g(f(oi)+{u, /(a))a) + (a, g(f(a) + (a, /(a))a))(a+ai+a 2 ) + (a, f(a))g(a) + 
(a, /(a)) 

(a + ai + a 2 ) = «i + gf(a) + (a,gf(a) + (a, /(a))g(a))(a + «i + 
"2) + (a, /(a)) (a + ai + a 2 ) = oli + gf(a) + (a, gf(a))(a + «i + a 2 ) = 
ai + a + 4>{A) + 0(B) + 90(a) + a + ai + a 2 + (a, 09(a)) (a + ai + 
a 2 ) + (a, d(j)(a))(a + «i + a 2 ) = a 2 + 0(5) + (a, 0(B))(a + aii + a 2 ) + 
0(A) + (a, 0(A)) (a + «i + a 2 ) + 90(a) + (a, 90(a)) (a + ai + a 2 ) = 
ftd'iax) + 9'0'(a 1 ). 

0'9'(a 2 ) + 9'0'(a 2 ) = 0'(e) + 0'(B) = a 2 + 0(B) + (a, 0(B)) (a + a x + 
a 2 ) + 0(B) + (a, 0(B)) (a + ai + a 2 ) = a 2 = a 2 + g'f'(a 2 ). 

<jf&{e) + 9'0'(e) = <j>'d{B) + 9'(a 2 + 0(B) + (a, 0(B)) (a + a x + a 2 )) = 
09(B) + d(a 2 ) + 90(B) = e + gf(B) = e + #7' (e) (recall that /(A) = 
f{B) = /(e) since /9 = in an AT-model). 

0'9'(a) + 9'0'(a) = <j/(d{<r) + (a,d(a))a) + (a, 9(<x))0'(a 1 + a 2 ) + 
9'(0(a) + (a, 0((r))a) + (a, 0((x))0'(ai + a 2 ) = 0(9(a) + (a, 9(<x))a) + 
(a,0(9(a) + (a, 9(cr))a))(a + ai + a 2 ) + (a, 9(cr))0(a) + 9(0(cr) + 
(a,0(a))a) + (a, 9(0(a) + (a, 0(a))a))(a + ai + a 2 ) + (a, 0(a))9(a) = 
°~ + gf{°~) + 0<9(°")) + ( a ? 90(a)))(a + ai + a 2 ) = a + gf(a) + (a, a + 
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gf(a)){a + a 1 +a 2 ). 



Besides, a + g'f'(a) = a + g'(f(a) + (a, f(a))a) + (a, f = 
o + gtf{o) + (a, f(cr))a) + (a, + («, /(o r ))a))(a + «i + a 2 ) + 

(a, f(a))(g(a) + a + a 1 + a 2 ) = cr + gf(a) + (a, gf(a))(a + ai + a 2 ). □ 

Observe that /i : C(if) — >■ C(H'), given by /i(a) = «i if a G if and 
/i(cr) = a for any a E H \ {a}, is a chain-group isomorphism. 



3. 3D Cubical Cohomology Ring 

In 12, 9[, the authors consider cubical complexes as the geometric build- 



ing blocks to compute the homology of digital images. In this section, we 
adapt to the cubical setting, the method developed in 0, for computing 
the simplicial cohomology ring of 3D binary- valued digital pictures. We must 



mention [14|, |10| as related works dealing with the cup product on cubical 
chain complexes in a theoretical context. 

Since we are working with objects embedding in M 3 , it is satisfied that 
homology and cohomology are isomorphic. However, cohomology has the ad- 
vantage over homology of having an additional ring structure given by the cup 
product, that is a topological invariant. This product provides information 
about the relationship between the generators of (co)homology, that enables 
to discriminate, for instance, pairs of cycles in different contexts, as in Fig. [TJ 
Notice that, in 3D, the only non-trivial cup products are those corresponding 
to elements of cohomology of dim. 1. If the cup product of two elements of 
cohomology of dim. 1 is not zero, then it is a sum of elements of cohomology 
of dim. 2. Recall that given an AT-model (P, H, f,g, <fi) for a polyhedral cell 
complex P, it is satisfied that H is isomorphic to the homology and to the 
cohomology of P. 

3.1. Cohomology Ring of Simplicial Complexes 

We recall now how the cup product is defined in the simplicial setting 
using AT-models: 

Definition 3.1. Q/ Let K be a simplicial complex. It is assumed that the 
vertices of K are ordered. Let (K,H, f,g,<p) be an AT-model for K. Let 



8 



. . . , /3 q } be the set of 2-simplices of H and let a\ and a 2 be two edges of 
H. The cup product of ct\ and a 2 is: 

g 

£((ai/ ^ <X2f)Wk)))Pk mod 2; 

where ct\f a 2 f on a 1-simplex (vi,Vj,Vk) with vertices Vi < Vj < 
is (ax, f(vi,Vj)) ■ (ct2, f(vj, Vk)); and (a±f ^ a 2 f) is extended to 2-chains 
(sums of 2-simplices) by linearity. 

Observe that for each k, g(/3k) is a sum of 2-simplices representing one cavity. 
Then, (a±f ^ Oi 2 f){g{l3k)) is a sum of Os and Is over Z/2 whose result is 
or 1. Therefore, J2k=i(( a if ^ a 2f){.g(fik)))fik is a sum of 2-simplices of 
H of dim. 2, representing the cavities obtained by "multiplying" the two 
representative cycles g(ai) and g{ct2) (think of the two tunnels of a hollow 
torus) . 

It is known that two objects with non-isomorphic cohomology rings, are 
not topologically equivalent (more precisely, they are not homotopic) [Til ] . 
To use the information of the cohomology ring for this aim, one can con- 
struct both matrices M and M' collecting the results of the cup product of 
cohomology classes of dim. 1 of each object, if the rank of M and M' are 
different, then we can assert that both objects are not homotopic (see 

3.2. Cohomology Ring of Cubical Complexes 

Now, let Q be a cubical complex. Our aim is to obtain a direct formula 
for the cup product on Q without making use of any triangulation. 

Suppose that the vertices of Q are labeled in a way that: 

(PI) Each square (vi,Vj,Vk,vg) of Q with vertices Vi < Vj < < ve has the 
edges (vi,Vj), (vi,Vk), (vj,vi) and (vk,vg) in its boundary. 

For example, a cubical complex whose set of vertices is a subset of Z 3 (the 
set of points with integer coordinates in 3D space 1R 3 ) with vertices labeled 
using the lexicographical order, satisfies PI. 

Definition 3.2. Let Q be a cubical complex satisfying PI and (Q, H, f, g, (p) 
an AT-model for Q. Let {f3i, . . . , (3 q } be the set of squares of H and let a>\ 
and «2 be two edges of H. The cup product of ct\ and a 2 is 

i 

^2((aif a 2 f)(g(Pk)))Pk mod 2; 
k=i 
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v k Vf V k V f V k V f V k V f V k tf ( 

Figure 4: Scheme of the cubical cup product. 

where a±f — 'q a 2 f on a square (vi,vj,Vk,ve) with vertices Vi < Vj < Vk < vt 
(see Fig. \4ty is: 

(«!, f(v t , vj)) ■ (a 2 , f(vj, v e )) + («!, f(v h v k )) ■ (a 2 , f(v k , v e )); 
and (cci/ ^ Q ^2/) extended to 2-chains (sums of square) by linearity. 

For simplicity, we sometimes use the notations (a ^ Q a ')(P) f° r («/ 
a' f)(g((3)), and analogously for the simplicial cup product. 

Example 3.1. Let Q be an abstract cubical representation of the hollow 
torus given in Fig. Consider the AT-model (Q, H, f, g,<p) for Q, given 
in Example \2.1\ Recall that H = {vo, (vo, v 2 ), (i>o, v±), (vo, v 2 , V4, v$)}; and 
g(v ) = vo, g(vo,v 2 ) = {v ,vi) + (vi,v 2 ) + (v ,v 2 ), givoiVi) = (v ,v 3 ) + 
(w3,f 4 ) + ( 1*0,^4) an d g( v o, v 2, v 8 ) is the sum of the squares of Q, repre- 
senting the connected component, the two tunnels and the cavity. 

Apply the formula given in Def. \3.S\ in order to obtain the cup product of 
(t> , v 2 ) and (v , w 4 ) in H : 

((v ,v 2 )*f ^ Q (v ,va)* f)(g(v ,v 2 ,v 4 ,v 8 )) 

■■= ((v ,v 2 ), f(v ,v 2 )) ■ ((v ,v 4 )J(v 2 ,v 8 )) 

+{(vo, v 2 ), f(v , v 4 )) ■ ((v , v 4 ), f(v A , u 8 )> 
= 1.1 + 0-0 = 1. 

then, (v ,v 2 ) ^ Q (1*0,^4) = (^0)^2,^4)^8)- Recall that (v ,v 2 ,v 4: ,v 8 ) is the 
square in H representing the cavity of the hollow torus. Therefore, the product 
of the two tunnels of the hollow torus is the cavity. 

The following theorem shows the validity of the definition of ^ Q (Def. 
13. 2p . That is, it is stated that we obtain the same result by applying the 
formula of Def. 13.21 to compute the cup product on the cubical complex, 
than making first a triangulation in order to obtain a simplicial complex, 
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Ml 101 




a) b) c) d) 

Figure 5: From a) to c), successive subdivisions; d) a cube subdivided in 6 tetrahedra. 

and applying the classical definition of the cup product given in Def. 13.11 
afterwards. 

Consider successive subdivisions of each cube of a given cubical complex 
Q until each one is converted in six tetrahedra, and such that each square 
(vi, Vj, Vk, Vi) of Q is subdivided by the edge (vi,ve) (see Fig. |5jd). Let 
us denote this resulting simplicial complex by Kq. Observe that with this 
particular subdivision, if (v p ,v q ,v r ) is a 2-simplex of Kq, with v p < v q < v r , 
obtained by a subdivision of a square of Q, then (v p ,v q ) and (v q ,v r ) will 
correspond to edges in Q. 

Theorem 3.1. Let (Q, H, f, g,4>) be an AT -model for Q. Let a and a' be 

two edges of H and G H a square. Let (Kq,H', f',g',(j)') be the AT-model 
for Kq obtained after successively applying Lemma \2. 21 Then, 

(a ^ Q a')(P) = (a ^ a>)(h(P)) 

where ' is the simplicial cup product given in Def. \3.1l ~-^ Q is the cubical 
cup product given in Def. \3.2\ and h : C(H) —> C(H') is the isomorphism 
defined at the end of Section [H 

ProofObserve that since a and a' are edges of H, then a, a' G H', that is, 
h{pi) = a and h{a!) = a'. Observe also that f(a) = f'(a) on any edge a since 

i /(«)• 

We have to prove that (af ^ Q a'f)(g(/3)) = (af — a' f')(g'(h(/3))). 

Let (3 = (vi,Vj,v k ,v e ). Let = (vi,Vj,Vt) and f3 2 = (v h v k ,v t ) be the two 
triangles obtained after subdividing f3 by the edge e = (vi,ve). Remember 
that h{0) = f3i and g'{h(f3)) = g'{Pi). Notice that coincides with g{0) 

if we replace (3 by f3± + j3% in the expression of g(/3). Therefore, it is enough 
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to prove that (af - Q a'f)(/3) = (af - a'f')(f3 1 + f3 2 ): 

(af ^ Q a'f)(/3) = (af — Q a'f)(vi, v d , v k , v e ) 
= (a, f(vi, Vj)} ■ (a 1 , f(vj, vt)) + (a, f(v h v k )) ■ (a', f(v k , v e )} 
= (a, f'(v h vj)) ■ (a 1 , f'(vj, v t )) + (a, f'{v h v k )) ■ (a', f'(v k , v £ )} 
= (af'^a'f)((3 1 + ^). 

This concludes the proof. □ 



4. Cubical Cohomology Ring of 3D Digital Pictures 

In this section, we develop the main bulk of the paper: beginning from 
a cubical complex Q that represents a 3D binary digital picture whose fore- 
ground has one connected component, first we compute an AT-model for the 
boundary dQ of the object; then, having in mind that the homology of dQ 
contains the homology of Q, we obtain an AT-model for Q with the rep- 
resentative cycles of homology generators lying in dQ; finally, applying the 
formula given in Section [31 the cohomology ring is computed from such an 
AT-model. 

4-1. From Digital Pictures to Cubical Complexes 

Each point of Z 3 can be identified with a unit cube (called voxel) centered 
at this point, with facets parallel to the coordinate planes. This gives us an 
intuitive and simple correspondence between points in Z 3 and voxels in R 3 . 

Consider a 3D binary digital picture I = (Z 3 , 26, 6, B), where B (the 
foreground) is finite, having Z 3 as the underlying grid and fixing the 26- 
adjacency for the points of B and the 6-adjacency for the points of Z 3 \ B 
(the background). We say that a voxel V is in the boundary of / if V G B 
(i.e., if the point of Z 3 identified with V is in B) and V has a 6-neighbor in 
Z 3 \R 

Take the cubical complex Q for B whose elements are the unit cubes 
(voxels) centered at the points of B together with all their faces. Observe 
that this cubical complex with vertices labeled by the corresponding cartesian 
coordinates and considering the lexicographical order, satisfies (PI) (see page 

ED- 

Without lack of generality, we consider that the foreground is connected. 
The elements of dQ are all the squares of Q which are shared by a voxel 
of B and a voxel of Z 3 \ B together with all their faces. 
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a) 



b) 



c) 



V 



Figure 6: a) A 3D binary digital picture I = (Z 3 , 26, 6, B); b) the cubical complex Q for 
B; c) the set of squares of dQ. 



4.2. AT-model for dQ 

Our interest now is to adapt the incremental algorithm for computing an 
AT-model given in 0, lEJ to the particular complex dQ. 

First, consider the set of edges and vertices of dQ as a graph and compute 
a spanning forest T. Let T 1; . . . , T m be the trees of T corresponding to the 
connected components of dQ. Fixed i, i = 1, . . . ,m, take a vertex Vi of T; 
and consider it as the root of Tj. 

Computing an AT-model (dQ, H, f, g, <fi) for dQ. 



INPUT: The complex dQ , 

the set {Ti, . . . ,T m } of trees of a spanning forest T of dQ , 
the set {v\,...,v m } of roots of the trees of T. 

1. Initialize f(a) := a, 0(a) := for any a G dQ ; H := {v±, . . . , v m } , 
U := {v : v is a vertex of dQ}, 

f(v) := Vi if v is a vertex of Tj for some i, i — 1, . . . ,m. 
For i = 1 to m do 

From £ = 1 to the height of Tj do 

For each vertex v at level £, and edge a linking v 
with its parent w do 

:= a + 0,H, U:=UU{a}, f(a):=0. 

2. While there are edges in dQ \ U do 

If there is a square c G <9Q \ U with exactly one edge 
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a G dQ \ U in its boundary do U :— U U {c, a} , 
/(a) :=/(0(c) + a), 0(a) := c + 0(<9(c) + a) , /(c) :=0. 
Else take an edge aedQ\U then # := if U {a} , U := U U {a} . 

3. While there is a square c in dQ\U do [/ := U U {c} . 

If /9(c) = then H := H U {c} . 
Else take an edge a in /(9(c) then H:=H\{a}. 
For each edge 6 in dQ \ T do f(b) := f(b) + (aj(b))fd(c) , 
0(6) := 0(6) + (a,/(i))(c + #)), /(c) = 0. 

4. For each er G H do o(cr) := cr + <f)d(a) . 
Output: the AT-model (dQ,H,f,g,<p) for 

The auxiliary set £/ is defined to indicate the cells which have already 
been used. In Step 1, neither the vertices nor the edges of Tj create cycles 
except for the root V{. In Step 3, if a square has edges in its boundary 
that created cycles in a previous step, then one of these cycles is destroyed. 
Otherwise, this square creates a new cycle (a cavity). In the last step, the 
representative cycles of homology generators are computed. 

Observe that all the steps of Alg. 14.21 are quadratic in the number of 
elements of dQ (worse case complexity) except for the last part of Step 3 
which is cubic in the number of edges of dQ \ T. 



Example 4.1. The AT-model (dQ, H, /, g, <fi) of a hollow cube dQ (see Fig. 
0) is: 





dQ 


f 





H 


9 


Step 1 


V 


V 





V 


V 




Vi, i = l,...,7 


V 


7(«i,v) 






Step 2 


Any edge b G Q\T 





Cb 






Step 3 











C 



where 7( Ui)i ,) is the only path in T from Vi to v; Cb is the square at which the 
arrow corresponding to an edge b in Fig. \7\c points, except for C(„ 4i „ 6 ) which 
is (i>4, v 5 , v e , Vj) + (vi, v 3 , v 5 , Vf); and C is the sum of the six squares of dQ. 
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a) b) c) d) 

Figure 7: a) A hollow cube dQ; b) a spanning tree T with root v; c) the "paths" q,; d) 
the cells of H. 



Face Reduction Process. 

INPUT: A cubical complex Q. Initially, K := Q . 

While there exist a, a' G Q \ dQ such that a 1 is in d(a) do 

For each cell c G K such that a' is in d{c) do 
redefine d(c) as <9(c + a) . 

Remove a and a' from the current K ; 
Output: the cell complex A". 



/ 






/ 


✓ 






/ 



1/ 



a) 



b) 



Figure 8: a) A cubical complex Q composed by two cubes c and <r sharing a square a' (in 
red) and all their faces; b) the squares of d(c) := d(c + a) in Q' . 
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4.3. AT-model for Q 

Now, we use a face reduction technique (see, for example, 0, 0, fli! ]) 
in order to obtain a cell complex K such that homology, cohomology and 
cohomology ring of K coincide with that of Q, and such that the cells of dQ 
are also cells of K. 

Observe that after the face reduction process, we obtain a cell complex K 
with the same topological information as Q but with much less cells. Now, 
starting from an AT model for dQ, compute an AT-model for K adding the 
cells of K \ dQ incrementally as follows: 

AT-model for K. 

INPUT: An AT-model for dQ: (dQ, H, f, g,(f>) and the cells 

{o~i, . . . , a m } of K\dQ ordered by increasing dimension. 
Initially, f K (a) := f(a), (f> K (a) := <j>{&) , g K (a) = g(a) for each a G dQ; 

f K (a):=0, <f) K (a):=0, for each aeK\dQ; H K := H . 
For i — 1 to i — m do 

take a cell, a, of f K d(o~i) , then 
H K :=H\{a}, " 
For k = 1 to k = i — 1 do 

UM ■= fxM + i^,f K {a k ))f K d{ai) 
4> K (o- k ) := (f) K (a k ) + (a, f K ((J k )){oi + 0k<9(>;)) 
Output: the AT-model (K, H K , f K , g K ,(j) K ) for K. 



Observe that in the algorithm, a cycle is never created because the cycles 
of dQ are also cycles of K. Therefore, when a cell U{ of K is added, then 
j K d[pi) is never null and therefore a class of homology is always eliminated 
(that is, a cell a of f K d(o~i) is removed from H K ). Alg. 14.31 is C(m 3 ), where 
m is the number of cells of K \ dQ (worst case complexity). 

4-4- From Cubical Complexes to Digital Pictures 

Now, given a digital picture /, suppose that we have computed an AT- 
model (K,H K ,f K , 

9k-, 0k) f° r K, following the steps given in Subsections |4JJ 
14.21 and H~3l For each a in H K , g K (a) is a representative cycle of a homology 
generator of K and, therefore, of Q, since the homology of K and Q coincide 
and the representative cycles are in dQ. Recall that if a is a vertex, then 
g K {o~) is a vertex representing a connected component; if a is an edge, then 
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a) b) c) 

Figure 9: a) A cubical complex dQ and a representative cycle of a homology generator; b) 
in green, voxels considered the first time the edges of the cycle are visited. In red, voxels 
considered in the second visit; c) in black, voxels representing the cycle. 




a) b) 
Figure 10: a) Two non-linked circles; b) two once-linked circles. 

g K (a) is a sum of edges representing a tunnel; and if a is a square, then g K (a) 
is a sum of squares representing a cavity. 

Given a representative cycle g K {c) of a homology generator, our aim in 
this subsection is to draw the equivalent cycle in the picture /. 

First, if g K {(j) is a vertex then, g K {&) is a face of a square in dQ. This 
square is shared by a voxel V of B and a voxel of Z 3 \ B. Then, associate 
the voxel V to the vertex g K (a). 

Second, if g K (a) is a sum of edges, suppose that g K (a) is a simple cycle 
(if not, it can always be decomposed in simple ones). Visit all the edges 
of the cycle in order. If an edge, a, and the next edge, b, are facets of a 
square o £ dQ, then associate the single voxel V of B which has a in its 
boundary, to the edges a and b. After that, visit all the edges that have not 
been associated to any voxel. If a voxel V is associated with the next edge of 
the current one, a, and there is a voxel V' £ B having a in its boundary, such 
that V is in the boundary of / and V and V are 6-neighbor, then associate 
V to a. If not, look at the previous edge and do the same procedure. 
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Q't 


Oil W Q Oil Oil W Q «2 «2 Ot 2 


Pi 

02 






Q'2 


a'l — Q a[ a\ ^ Q a' 2 a' 2 — Q a' 2 


P'l 


1 
1 



Figure 11: In yellow, representative cycles of the two tunnels of Q[; in green, representative 
cycles of the two tunnels of Q' 2 ; at the bottom, the tables for the cup product on Q[ and 
Q' 2 , where a% (resp. aj), i = 1,2, are representative cycles of the two tunnels and Pi (resp. 
/?,'), £ = 1,2, are representative cycles of the two cavities of Q[ (resp. Q' 2 )- 

If not, take any voxel of B that contains a, having a 6-neighbor voxel in 
Z 3 \B. 

Finally, if g K (a) is a sum of squares, associate, to each square a, the single 
voxel V of B which has a in its boundary 

4-5. Cohomology ring of Q 

Given a digital picture (Z 3 , 26, 6, B), its associated cubical complex Q and 
having computed an AT- model (K, H K , /k^kj^k) f° r K, the last step of the 
process is the computation of the cohomology ring. This can be performed 
using the formula for the cubical cup product given in Def. 13.21 

Example 4.2. This example shows an application of — y q to discriminate 
different embeddings of the same object. Consider the cubical complex Qi 
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Figure 12: A configuration of 6 linked circles and the results of the computation of the 
cup product of the cubical complex associated to the white voxels of the picture. 

(resp. Q2) associated to a digital picture where the set B consists in two 
once-linked "circles" (resp. two unlinked "circles"). See Fig. [73 Both 
complexes have two tunnels and no cavities, so these properties are not able 
to distinguish them. 

Now, denote by Q[ and Q' 2 , the cubical complexes associated to the back- 
ground of Ii and I2 (white voxels of Fig. QJJj). Compute an AT -model for Q\, 
and its cohmology ring, fori = 1,2. We obtain that the multiplication table 
for the cup product on Q[ is null whereas on Q' 2 is not (see Fig. [77]) . This 
fact allows us to assert that the two complexes Q[ and Q' 2 are not topologically 
equivalent. 

Example 4.3. Consider the picture in Fig. [J_H The cubical complex associ- 
ated to the white voxels of the picture has 1 connected component, 6 tunnels 
and 3 cavities. The cup product is trivial for any two pairs of homology 
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al pha 1 



al pha 2 



al pha 3 



al pha 4 
al pha 5 



al pha 6 



al pha 



alpha 1 alpha 2 alpha B alpha 4 alpha 5 alpha 6 alpha 7 



Figure 13: A configuration of 7 linked circles and the results of the computation of the 
cup product of the cubical complex associated to the white voxels of the picture. 



classes except for two tunnels named as 012 and wich is the sum of two of 
the three cavities (see Fig. [JDj] . 

Finally, Consider the picture in Fig. [73]. The cubical complex associated 
to the white voxels of the picture has 1 connected component, 7 tunnels and 
12 cavities. The results of the computation of the cup product can be seen in 
the table of Fig. [73J where "CP i j " means the sum of the cavities i and j . 



Example 4.4. MRA of chest showing the heart and great vessels. 

Originally it was a stack of 14 gray-scale digital images in DICOM format, 
each one of size 320x320. This 2D stack can be see as the 3D digital image 
"chest" of dimensions 320x320x14 (see Fig. 14 a). After a binary process, 
we obtain the picture I in Fig\T^b). "resized" 80x80x28 

Let Q be the cubical complex that represents I , let dQ be 
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Number of black voxels of I 18, 062 

Number of cells of Q 179, 769 

Number of cells of dQ 66, 143 

Number of cell of K 

AT-model implementation execution time 



5. Conclusions and Future Work 

In this paper we present formulas to directly compute the cohomology 
ring of 3D cubical complexes and develop a method for the computation on 
3D binary-valued pictures. This computation on cubical complexes can be 
regarded as a starting point to compute the cup product on general polyhe- 
dral cell complexes, which is, in fact, our last goal. As related work, we must 
mention [gj], where homology of 3D pictures using particular cell structures 
provided by the 26-adjacency is performed. The restriction to the 3D-world 
allows to work over Z/2, what facilitates the calculus. However, a harder 
task could be the one of extending the formulas of the cohomology ring to 
higher dimensions what could be applied to more general contexts out of 
digital images. In this sense, cohomology ring of nD simplicial complexes 
using AM-models and working in the integer domain, has been established 
in [6|. Another goal for future work is the one of applying theoretical results 
to irregular graph pyramids and compute the cohomology ring on the cell 
complexes associated to such structures. In the paper 0], representative co- 
cycles for cohomology generators on irregular graph pyramids are computed, 
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what can be considered a first step in this direction. 
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